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Abstract. This paper describes the apphcation of a laser diffraction technique to the study of electrocon- 
vection in nematic hquid crystal cells. It allows a real-time quantitative access to pattern wave lengths and 
amplitudes. The diffraction profile of the spatial periodic pattern is calculated and compared quantita- 
tively to experimental intensity profiles. For small director tilt amplitudes the phase grating generated 
in normally incident undeflected light and the first order term correction from light deflection is derived 
analytically. It yields an / oc (/?^ dependence of the diffracted intensity / on the amplitude of director de- 
flections. For larger director tilt amplitudes, phase and amplitude modulations of deflection of light in the 
inhomogeneous director fleld are calculated numerically. We apply the calculations to the determination 
of the director deflection and measure growth and decay rates of the dissipative patterns under periodic 
excitation. Real time analysis of pattern amplitudes under stochastic excitation is demonstrated. 

PACS. 4 2.70.Df (Liquid Crystals), 
47.20.-k (Hydrodynamic instability), 

78.20.-e (Optical Properties of bulk materials and thin films). 



1 Introduction 

" Electrohydrodynamic convection (EHC) in nematic liq- 
uid crystals is one of the standard systems of dissipative 
pattern formation. It has been studied extensively during 
past decades. As a consequence of anisotropic properties 
of nematic phase, the system is particularly rich in mor- 
phology. Among the advantages of this system for exper- 
imental characterization are the easy control of electric 
excitation fields, convenient time scales and the straight- 
forward observation techniques. 

The equations describing the fundamental mechanism 
yield two dynamic regimes: conduction and dielectric struc- 

' tures. In addition to the primary instability toward sim- 
ple roll patterns with wave vectors normal or inclined to 
the preferential alignment of the director (optic axis) of 
the system, a variety of secondary instabilities have been 
described. Besides the investigation of arrays of parallel 
rolls, scientific interest recently focussed on defect struc- 
tures and localized convection states 1 2j^||l||^. The most 
successful and widely used method for the investigation 
of the complex convection patterns is optical microscopy, 
based on the shadowgraph method jHI- It utilizes the de- 
flection of light rays in the spatially modulated director 
field of the nematic. It has been applied to the determina- 
tion of wave vectors, onset thresholds and subcritical fluc- 
tuations of convection. One of the problems encountered 



in this very efficient observation technique is the complex- 
ity of the optics in the periodically deformed director field. 
Simulations of the optical profiles have been presented by 
several authors pi7ll51IMTnilllU12j . and the consequences 
of in-plane director twist have been considered 13 , 14 . Al- 
though light propagation in such a two-dimensionally in- 
homogeneous medium has been treated theoretically with 
different approximation methods, and the qualitative re- 
lation between director structure and observed intensity 
profile in the microscope is well established, the method 
fails to provide quantitative access to deflection ampli- 
tudes of the director field. With varying amplitudes of the 
spatially periodic director modulations, both the positions 
of focal planes of the patterns and the intensity profiles 
at given focal planes of the microscope change in a com- 
plex way. Thus, the power of the method lies primarily in 
a quantitative determination of the pattern wavelengths 
and the topology of defect structures. Moreover, a fast 
(real-time) observation of the pattern amplitude dynam- 
ics requires considerable bandwidth and signal processing 
speed. 

The wave vector spectrum can be determined quali- 
tatively from a Fourier transform of digitized microscopic 
transmission images 12,15,16 . A more efficient, quanti- 
tative way is the quasi-optical Fourier transformation by 
means of laser diffraction. Laser diffraction has been ap- 
plied earlier to the study of Williams domains or com- 
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parable dissipative patterns of nematics by Akahoshi et 
al. Vistin and Yakovenko and Miike et al. jl9j . 
In these studies, the evaluation of the scattering profile 
remained to a large extend qualitative. A first experimen- 
tal and theoretical study of the laser diffraction efficiency 
of EHC has been presented by Carroll [2^, and Kash- 
now [21 . Scattering spots designated to gratings in phase 
and amplitude of transmitted light can be distinguished. 
They are generated by spatially modulated optical path 
and light deflections, respectively. 

In a structurally similar system, Bouvier and Scharf 
have employed the Jones matrix method to calculate the 
diffraction efficiency of periodically deformed director fields 
in cells with structured electrodes and compared it to ex- 
perimental data [22j . Their method depends upon the as- 
sumption of normal undeflected light transmission through 
the medium and describes only the phase grating. 

Comprehensive analysis of diffraction gratings formed 
by EHC has been presented by Zenginoglou, Kosmopou- 
los and Papadopoulos j9l|23li24ll25il26l . Various aspects of 
laser diffraction by EHC have been considered, like the 
test of the validity of geometrical optics ^5 , diffraction 
under oblique incidence 23 , director oscillations and re- 
laxation, and the dielectric regime 26ll27j . In a study of 
stochastically excited EHC 28 ,29 ,30 , laser diffraction has 
been successfully applied to characterize fundamental scal- 
ing laws in the statistical description of pattern dynam- 
ics. The advantage of the laser diffraction technique over 
shadowgraph observations is particularly evident in such 
an experiment where data reduction is necessary to pro- 
cess pattern wavelengths and the trajectory of the pattern 
amplitude in real time. 

The previous study of fundamental scaling laws in stoch- 
astically driven EHC was based on several properties of 
the diffraction profiles which have not been explicitly given 
there 28,30 . This manuscript deals with the underlying 
optical principles, it provides a justification of the quan- 
titative relations between diffraction profiles and director 
field deflection amplitudes, and moreover, derives the ab- 
solute diffraction efficiencies which allow the determina- 
tion of not only growth rates and relative pattern ampli- 
tudes but also the director tilt amplitudes in a quantitative 
way. We demonstrate the application of the laser diffrac- 
tion technique to the study of amplitude dynamics of de- 
terministically and stochastically excited patterns. A rig- 
orous treatment of light propagation in two-dimensionally 
inhomogeneous director fields of EHC has not been pre- 
sented so far, therefore we will discuss the validity of sev- 
eral approximations. We recollect the methods to calculate 
light propagation and the corresponding phase and am- 
plitude modulations in two-dimensionally inhomogeneous 
director fields. A weakly nonlinear analytical calculation is 
compared with numerical simulations of the full nonlinear 
equations and with quantitative experimental data. 

The paper is organized as follows: a short introduc- 
tion into the basic principles of the Carr-Helfrich mecha- 
nism and the involved dynamic equations is given in the 
second section. We introduce the experimental setup and 
the qualitative structure of the diffraction patterns in the 



third section. In Sec. 4, we derive the analytical formula for 
the phase and amplitude gratings generated by a weakly 
distorted director field, which also accounts for first order 
effects of light deflection. In the course of the fourth sec- 
tion, the quantitative numerical calculation of the diffrac- 
tion profile is performed. Although a great part of the 
equations derived in this section have been communicated 
in earlier work by other authors, we consider it helpful 
to include a comprehensive treatment of the optical back- 
ground here, in particular because in literature sometimes 
there seem to be contradictory details of the calculations 
(see below). One obtains numerically the phase and am- 
plitude gratings produced by the mesogen layer for nor- 
mal and oblique incidence of monochromatic light. This 
allows quantitative predictions from the combined effects 
of ray deflection and optical path length modulation. We 
describe the diffraction efficiency of periodic nematic di- 
rector structures and compare our calculations with the 
approaches proposed in literature. The numerical and an- 
alytical results are tested by comparison with the experi- 
ment. Finally, we apply the method to the determination 
of growth rates and Lyapunov exponents of pattern ampli- 
tudes in EHC and demonstrate the power of the method 
to determine real-time amplitude fluctuations of director 
field modes. 



2 Electroconvection 

Electrically driven convection in liquid crystals bases on 
the interaction of free charges in the mesogen with exter- 
nal electric fields, and the coupling of fluid flow to the 
deflection of the nematic director. A comprehensive re- 
view is given, e.g., in Refs. |3Tll32ll33ll3^ . The essential 
variables describing the structures are the spatially mod- 
ulated charge distribution q{x^ z) and the director tilt 
angle (^(x, z), both are coupled via the electrohydrody- 
namic equations. In an oscillating excitation field, the two 
quantities have qualitatively different dynamic behaviour. 
For the diffraction experiment as well as for the conven- 
tional shadowgraph images, only the director field modu- 
lation is relevant and accessible. 

A sketch of the experimental geometry is given in Fig.^ 
The nematic director alignment at the glass plates is fixed 
along X by surface treatment. The ground state is a uni- 
form director field in the cell. When an electrical field E = 
U/d is applied between the transparent ITO-electrodes at 
the glass plates, free charge carriers (ionic impurities or 
dopants) in the nematic fluid are accelerated and initiate a 
macroscopic flow. The conductivity anisotropy of the ma- 
terial in combination with small fluctuation modes of the 
director tilt lead to lateral charge separation in x direction 
and a periodically modulated flow field, which in turn 
couples to the director field by hydrodynamic equations. 
At the critical field Ec^ stabilizing elastic and dielectric 
torques on the director are outmatched by destabilizing 
hydrodynamic torques. Under standard conditions, the 
system exhibits a forward bifurcation to normal rolls (wave 
vector along the easy axis of the director) or to oblique 
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Fig. 1. Schematic drawing of convection rolls and director 
field in a nematic sandwich cell. A snapshot of the spatial mod- 
ulations of director and charge fields ((^, q) in the cell midplane 
is sketched. 



rolls. Threshold voltage and critical wave number are fre- 
quency dependent. The pattern stability diagram of the 
two sandwich cells studied here is shown in Fig. [21 The 
nematic material is Mischung 5, a mixture of four di- 
substituted phenyl-benzoates 28 , material parameters in 
Table I. The first cell has been prepared with the pure, 
undoped material, which has a low conductivity and cor- 
respondingly low cut-off. 

The nematic mixture in sample 2 has been doped with 
0.5 mass%o tetrabutyl-ammonium bromide. Therefore its 
conductivities are much higher and the increased amount 
of charge carriers shifts the cut-off frequency out of the fre- 
quency range of measurements. It leads to a much more 
stable pattern amplitude characteristics near the thresh- 
old. In the undoped material, the content of charge carriers 
is comparably small ('natural' impurities after synthesis) 
and in the experiment, the threshold voltage is subject to 
certain small but measurable long-term fluctuations. 

In Fig. El one distinguishes the low frequency 'conduc- 
tion' regime and the high frequency 'dielectric' regime of 
cell 1, with a distinct jump in the wave number at the 
cut off frequency. Only the low frequency regime is in the 
accessible frequency range in cell 2. All measurements in 
this study are performed in the conduction regime, where 
the director field performs only moderate oscillations syn- 
chronous with the excitation frequency, but keeps its sign 
during the field cycles. However, there is no principal lim- 
itation for an application of the presented setup to struc- 
tures of the dielectric regime j ^l2 71135 j . 

The mathematical description is based on the Maxwell 
and Navier Stokes equations. The linear stability analysis 
of the torque balance uses a test mode ansatz for director 
deflections and charge density modulations 



(f{x^ t) = (ft cos{kxx) cos{kzz), 
q{x, z, t) = qt ^\n{kxx) cos{kzz), 



(1) 

(2) 



where kx is the periodicity of the pattern. Because the 
director is fixed parallel at the glass plates, the boundary 
conditions (p{^z = d/2) = enforce kz = (2n + l)7r/(i, 
n integer. Near onset, we consider only the ground mode 
kz = Ti/d. In the investigated parameter regions here only 
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Fig. 2. Stability diagram of the electroconvection patterns at 
the first instability measured for cell 1 (25.8 pm thick, ■) and 
cell 2 (48.5 pm thick, o) under square wave excitation. The 
thresholds under sine wave excitation are not much different 
from the square wave case (at comparable effective voltages). 
Lines are from the analytical calculation, where known mate- 
rial parameter from independent experiments, have been used, 
if available. The remaining unknown parameter have been ob- 
tained from the fit to the experimental data. The cut off fre- 
quency z/c = 51 Hz for cell 1 separates conduction and dielectric 
regimes. In case of cell 2 the much higher cut-off frequency is 
outside of the presented range. 



normal rolls appear, and no y dependence has to be con- 
sidered. Linearization leads to a linear ordinary differen- 
tial equation system in (g, ip). The solution for a constant 
electric field amplitude E involves a 2 x 2 non-symmetric 
time evolution matrix T jHEl which depends on the wave 
vector kx of the particular test mode 



{t) = T{E,kxMlm 



(3) 



At square wave excitation, where only the sign of E alter- 
nates, the time evolution at points with alternating sign 
of E is given by a product of matrices with constant co- 
efficients 



(t„) = T±(Z\t„) ■ ■ • T-{Ah)T+{Ah) (^Jj (0), (4) 



(5) 



where Ati are the time intervals between consecutive jumps. 
The solution at intermediate times is calculated with Eq. 

0. The involved material parameters are listed in Table 

1, further details are given in [36] , 

The largest of the two real eigenvalues of the matrix 
product in Eq. Q is related to the Lyapunov exponent of 
the system 36 and describes the asymptotic behaviour of 
a small initial perturbation in (g, Lp). In case of a periodic 
excitation, all Ati = At = l/{2v) are equal and after n 
full periods of the ac excitation, the product in Eq. 
can be split in repeated blocks of the product T+T~ 



(tn) = (T+T- 



w^g(o). 



(6) 
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For tn ^ At^ the amphtudes of both variables grow or 
decay exponentially, and the (dimensionless) largest real 
eigenvalue Xi{E^ kx^ i') of T+T~ gives the growth or decay 
rate vXi. The maximum \i{E^kx-,y) of all kx selects the 
critical wave number k^. The theoretical threshold field 
Ec at a given frequency v is determined by the first pos- 
itive value of Xi^E^kc^f) with increasing E. This value 
coincides with results from the Floquet theorem. 

With any optical detection methods (shadow graph or 
diffraction), only (^t is observable and asymptotically for 
At^ (in the limit of small ipt) 



vXi t 



(7) 



where, (po is the initial amplitude of the considered mode, 
related to fluctuations of director and charge fields. For 
large amplitudes, (ft is limited by nonlinear it ies that will 
not be considered in the linear treatment. For sample 1, 
the theoretical threshold curves Uc{i^) = E^d and /Cc(^) 
can be fitted to the experimental data over a wide fre- 
quency range with most of the involved material param- 
eters taken from independent experiments, see Table I. 
Some remaining unknown viscoelastic parameters are ob- 
tained from the fit. With this completed set of parameters, 
the theoretical Ai(£^, /cc, dependence can be calculated 
analytically as a function of the excitation frequency v. 



3 Laser difFraction experiment 

The diffraction of laser light by the director pattern pro- 
vides the opportunity to analyse the spatial mode spec- 
trum in real time. Figure 01 sketches the experimental 
setup consisting of a low power {^1 mW) He-Ne laser, 
the liquid crystal (LC) cell mounted in a Linkam micro- 
scope hot stage TMS 600, a photodiode or alternatively a 
diffusely reflecting screen for 2D camera images. 

The photodiode can be moved by a stepper motor in 
horizontal x direction across the scattering image. Its aper- 
ture is 3 mmx3 mm. At a distance I of approx 800 mm 
from the LC cell, this corresponds to an angular resolution 
of 5 mrad. The 2D images of the CCD camera are used 
for the qualitative characterization of the diffraction pat- 
terns only (see appendix), while all quantitative intensity 
measurements are performed with the photodiode. 

The test modes of Eqs. (|1I2|) correspond to a one- 
dimensional stripe pattern along x in the microscope im- 
age. For the sample cell 1 studied here, it appears at in- 
termediate frequencies, from the Lifshiz-point (below 20 
Hz) to the cut-off frequency 50 Hz), as the first insta- 
bility. At lower frequencies, the wave vector at onset has 
a non-zero y component. 

Figure ^ demonstrates the two cases of normal and 
oblique rolls for square wave excitation at 30 Hz and 10 
Hz, resp., at voltages slightly above the convection thresh- 
old. In such cases where one or two superimposed wave 
vectors form the spatial structure, the diffraction experi- 
ment provides the pattern wave lengths, orientations and 
amplitudes. Defects and domain sizes will primarily in- 
fluence width and fine structure of individual diffraction 
peaks but are not directly accessible from the profiles. 
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Fig. 3. Sketch of the experimental setup. The driving voltages 
are synthesized by a computer with digital analog (DA) con- 
verter and analog amplifier. The photodiode signal is sampled 
by the computer at a maximum rate of 7 kHz with an accuracy 
of 12 bit, or alternatively sampled by a digital voltmeter at a 
rate of 50 Hz with an accuracy of at least 6 digits. 
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Fig. 4. Snapshots of diffraction images (a,b) and respective 
microscope images (c,d) for normal rolls at 30 Hz (a,c) and 
oblique rolls at 10 Hz (b,d), cell 1. Numbers in (a) mark the 
diffraction order. In the lower left part in (c), a localized defect 
of the roll pattern is visible. The cell thickness is 25.8 pm. 
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4 Optics 

4.1 Light propagation 

The Ught intensity at the position of the diffraction re- 
flexes is directly related to the amplitude of the spatial 
director mode. The problem of light propagation in a pe- 
riodically modulated director field has been investigated 
by different groups before |Hlinil2 0112111^ . These previous 
works can be grouped into two different approaches. One 
is based on the solution of the Maxwell equations, the cal- 
culation of the spatial distribution of E and D. Yet, a 
complete solution of the Maxwell equation with bound- 
ary conditions can only be obtained numerically for the 
present problem. A linearized wave propagation approach 
can be found, for example, in Ref. 25 . The alternative 
method is based on the calculation of light ray's using 
the eikonal method jSZI oi" crystal optical methods 
Here, we give a short outline of the calculation of phase 
and amplitude of the laser light passing the LC layer using 
crystal optics. 

The polarization vector E of the incident laser light is 
adjusted along the director easy axis, since only the ex- 
traordinary wave is relevant for the diffraction effect (see 
Fig. Ordinary rays pass the LC layer without deflec- 
tions and phase modulation. 

In the lowest order of approximation, we may assume 
a straight propagation of the electromagnetic wave with- 
out any deflection of the light beam, r{z;xo) being the 
x-coordinate of the ray 

r {z]Xo) = dr{z;xo)/dz = (8) 
r{z;xo)=xo. (9) 

It yields no amplitude modulation but a first estimation of 
the phase profile of the light that has penetrated the cell. 
Under these assumption the phase tps{x) in dependence 
on the amplitude of director deflections is given by (see 
Fig.ISl) 

d/2 

2ps{x) = /cl J nl^{r, z) dz = Kn^d + /AV^s(^), (10) 

-d/2 

ny/3) = ""^'^ =, (11) 

J nl cos^ + si^^ P 

/3{r^ r', z) = (ft cos{kxr) cos{kzz), (12) 

2 2 

Z\V^s(^o) ^ -Kued-^^-^ ipl cos(2A:^xo), (13) 

where n^^ is the effective refractive index and no,ne are 
the ordinary and the extraordinary refractive indices of 
the nematic material (see Table 1). Only the first non con- 
stant term in the series expansion of n^^ enters the result, 
the constant phase does not contribute to the diffraction 
profile. 

The exact calculation of light propagation basing on 
the calculation of ray paths {r{z] xq), z) uses the Fermat 
principle. From the symmetry of the problem, no change 




Fig. 5. Definition of angles on the ray ellipsoid, Lp : director 
deflection, n /2 — (3 : angle between optic axis and Poynting- 
vector S, 7r/2 — : angle between optic axis and normal of 
plane wavefront k. 

of polarization of the light can occur. If the wave number 
of the laser light /cl is much larger than the wave number 
of the pattern, the Fermat principle can be applied in 
the birefringent material. The minimum condition is 

d/2 

j nl^(r, r' , z) ds = min, (14) 

-d/2 

n\^{r, r' , z) = \J cos^ /3 ^ sin^ /3, (15) 

/3(r, r\ z) = cpt cos{kxr) cos{kzz) — arctanr', (16) 

ds = vTTr^dz, (17) 

where nl^ is the effective ray index ^. Applying the Euler- 
Lagrange formalism leads to a ordinary differential equa- 
tion of second order in the displacement r{z]Xo) of a ray 

2 2 
T) — T) 

/' = (^,^^(l + rO(ti+t2), (18) 

ti = kz cos{kxr) sin{kzz) (n^ cos^ (3 — n^r' cos^ (3 sin (5— 
-nlr' cos (3 sin^ P - ^i^"^ P) ^ 

t2 = kx sm{kxr) cos{kzz) (nlr cos^ P -\- cos^ P sin /3+ 
-\-nl cos /3 sin^ /3 — n^r' sin^ 0) . 

Eq. (|T8j) is an exact result which we use in a Runge- 
KuTTA algorithm for a numerical computation of the ray 
paths. Expanding of (|14|) to first order in the amplitude 
(ft of the spatial mode of director deflection yields 

2 2 
Tl — Tl 

r" ^ -^^r^^tkz cos{kxr) sm{kzz). (19) 

The integration of Eq. (^HJ with the initial values 

r\z=-d/2 = xo, r'\z=-d/2 = (20) 

^ In literature, e.g. |()ll38ll39) . there has sometimes been a 
confusion about the appropriate usage of the ray index or re- 
fractive index in the computation of light propagation through 
a birefringent nematic layer with periodically deformed direc- 
tor field. 
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and the approximation cos{kxr) = cos{kxXo) leads to a 
path which enters at xq and 

1 — 77/^ 

r{z]Xo) =xo- ^ "" (ptil^ sin k^z) cos kxXQ. (21) 

kz 

To describe the focussing effect in the shadowgraph method 
the second order term in (ft must be considered, it can be 
found in Ref. 6 . For the diffraction pattern, mainly the 
phase is important (at least for small director deflections) 
and therefore no consideration of higher order terms are 
necessary. In order to calculate the resulting phase beyond 
the cell, we have to distinguish between the direction of 
propagation of energy flux S and the normal of wavefronts 
k (Fig.EJ. In uniaxial birefringent material the relation is 

m 

The phase can be calculated using different indices 

|- =1 nl^s = I <ffd|k| = I <gCos(/?-^)ds. (23) 

Finally, the calculated lateral phase difference along curved 
paths at the exit position for small director deflections is 
basically the same as the result for straight transmission 
(P|) 

A^{r,z)\,=d/2 = ^M^o)- (24) 

This result is in agreement with the relation in (371 calcu- 
lated from the eikonal equation. Also, the next order term 
in if^ is stated in ^j. The consideration of the displace- 
ment r{z;xo) — xo gives 

A^l){x, z)\^=d/2 - -Kried ^^^^ ^ cos(2kxx), (25) 
X = r(c^/2;xo), (26) 

X denotes the exit position r{d/2;xo) of a beam entering 
at Xq which may differ from xq due to the ray deflection. 
Therefore, the first order bend of the light paths is an ef- 
fect of the birefringence (the inclination of the optic axis 
of the nematic material respective to the incident beam 
direction) and is not the result of the periodic modula- 
tion of the refraction (connected with the ray index) in 
x-direction. 

The displacement of the rays after propagating the cell 
and the conservation of energy gives the intensity X in 
terms of the incident intensity 2o, and the amplitude of 
the electric field £ (with f ^ oc X) in the exit plane 

XqSxq = XSx^ (27) 
Xodxo =X [r{z;xo + dxo) - ^(^; ^o)]L=rf/2 ^ (28) 

X S^x) / k^nl-nl . , .o.x 
- = ^=(l + 2---^^.sm/...j . (29) 

The ansatz (|T9|) is only correct in the limit of small (ft 
when rays leaving the cell in parallel in good approxima- 
tion. A correction resulting from more exact treatment 
based on VS = can be found in Ref. 




-12,5 0,0 37,5 

x-ax\s [[im] 



Fig. 6. Numerically (solid) and analytically (dotted) calcu- 
lated ray propagation (a), intensity (b), and relative phase (c) 
of light after penetrating the nematic layer at 2; = 25 pm. The 
assumed director tilt {cp = 20°) is schematically depicted in 
(a), the grey scale visualizes the effective refraction index for 
straight light propagation. The periodicity of the amplitude 
grating is that of the director field. In contrast, the dominating 
phase modulation has twice the wave number of the director 
field. A cell thickness of 50 pm and optical parameters of Tab. 
I have been assumed. 

Figure EK) illustrates the light propagation from dif- 
ferent entry positions along the x-axis in the case of a 
strong director tilt amplitude (p = 20°. It compares the 
numerically calculated paths with the analytical result 
from (|21|) . Also the intensity and the phase at exit po- 
sition are depicted in Fig|Sl3,c). The analytical result for 
the phase differs hardly from the exact numerically calcu- 
lation whereas the first order approximation in intensity 
differs even qualitatively from the numerical values. Fortu- 
nately, the diffraction efficiency is dominated by the phase 
profile and therefore the conclusions drawn from diffrac- 
tion profiles about director deflections are correct up to 
large deflections angles. 



4.2 Diffraction profiles 

The cell is illuminated with a normally incident planar 
wave, polarized in x direction. The area contributing to 
diffraction is a circular spot with radius s ~ 0.5 mm. In 
general, the electric field at the rear of the cell can be 
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(30) 



written as ^ 

£\z=d/2 = £{x)e'- 

with the amplitude E{x) and AiIj{x) the lateral phase dif- 
ference of the wave at the position (x, z = d/2). For a one- 
dimensional modulation (Eqs. (|1I2|) ). each location x can 
be considered as the origin of a spherical wave. The diffrac- 
tion intensity at ^ ^ 5 into the angle in the (x, z) plane 
is 



e 

hi- 



xki^ sin 9, 



(31) 
(32) 



where /cl is the wave number of the incident light, the 
vector 1 connects the cell with the detector position and ^ 
the distance between cell and detector. Integration over a 
circular area with radius s of the illumination spot gives 
the amplitude of the complex wave with the corresponding 
intensity 



s \j8^-y'^ 

E{0)(x J J f(x)e^[^^^'^"^+^^(^)]dxdy, (33) 

m^\E{e)\\ (34) 

To calculate the complete diffraction function numer- 
ically, we consider both the spatial modulation of £{x) 
(amplitude grating) as well as the phase modulation AiIj{x) 
(phase grating). Basically, the first one is effective for 
shadow graph images. In contrast, when the director mod- 
ulation is small, the latter plays the dominating role in the 
diffraction characteristics. Therefore, we will consider for 
analytical approximations the phase grating alone. 

In case of small director deflections, the differences be- 
tween exit position x and entry point xq can be neglected 
and the periodicity of the phase grating can be written as 



Atlj{x) = Z\V^maxCOs(2/c^x). 



(35) 



It is twice that of the director field and therefor diffrac- 
tion reflexes from the the phase grating appear only at 
even order n. On the other hand the periodicity of the 
amplitude grating (|29|) is the same of the director field 
and it contributes also on odd order reflexes. 

The integration of ((33j) gives the intensities In = I{On) 
of the n— th order diffraction spots. With the assumptions 
£{x) =const. and (pTHj) the intensities at angles 6n are de- 
scribed by Bessel functions Jn/2 with the amplitude of the 
laser light phase modulation in the argument 







2^[(n/2)! 



where /q is the intensity of the mean beam at the ground 
state. Using (|24|) and (|13|) gives the final quantitative rela- 
tion between normalized diffraction intensity and director 
tilt amplitude at the dominating second order reflex 



8n2 



(f^, where h <. lo- (37) 



The numerically obtained diffraction efficiencies in the 
limit of small cpt confirm the relation (|37|) including the 
prefactor. The analytical approximations are satisfactory 
up to (p ~30° for thin cells. 

The numerical calculation offers an easy way to con- 
sider effects from oblique illuminations, an important point 
to understand the sensitivity of the diffraction images to 
a non perfect sample orientation. A detailed discussion is 
found in the appendix. 



4.3 Experimental test 

We tested the presented calculations with the setup sketch- 
ed in Fig.jSlwith sample 2. The temperature is stabilized 
at 32 °C and the measurement is performed at 500 Hz sine 
excitation. Due the high conductivity the cut-off frequency 
is shifted above 600 Hz and the pattern is more perma- 
nent at 500 Hz sinodial excitation than in cell 1. A circular 
aperture with 0.5 mm in diameter defines the illumination 
spot. A smaller movable photo diode (0.3 x 0.6 mm) is po- 
sitioned in a distance of ^ = 820 mm, to obtain a high res- 
olution in diffraction angle. The analog-digital-converter 
is a programmable Kethley multimeter with a resolution 
of 24 bit and a sampling rate of 50 s~^. 

Figured shows the baseline corrected intensity along 
X direction in the x-z plane {y = in Fig. ^ for two 
applied voltages slightly above the instability threshold. 
All profiles are normalized with the transmitted primary 
beam intensity at ^ = . Note the logarithmic intensity 
scale. In order to discriminate the diffracted light from 
a small constant background (scattered light from glass 
plates and small amplifier offset), we subtract the constant 
signal of the order of 10~^, detected at large deflection an- 
gles {6 > 0.5 rad), from all measurements. In superposition 
with the constant offset, small fluctuations in order of 10~^ 
are observed which drop below the detection level when 
the sample temperature is increased above 80 °C into the 
isotropic phase. The intensity profiles are compared with 
numerical calculations for different director amplitudes. 
The best fits lead to director deflections (f = (7.4 ± 0.3)° 
and if = (12.6 ± 0.3)°, e = U/Uc - 1 2 x 10^ and 
e ^ 5 X 10^, resp. (Fig.EJ. The normalized intensity at 
the second order reflex coincides well with the value ob- 
tained from analytical treatment (Eq. (pTZ)) ). 



5 Applications 



(36) 5.1 Study of EHC near the threshold 



The sensitivity in the detection of small pattern ampli- 
tudes and the quantitative relation between diffraction in- 
tensity and director tilt provides the opportunity to study 
electroconvection near the onset threshold experimentally. 
Amplification of thermal fluctuations slightly below the 
threshold has been studied previously with the shadow- 
graph method 



8 



Th. John et al.: Laser diffraction by EHC 




Time t [s] 

2000 



-0.09 -0.06 -0.03 0.00 0.03 0.06 

diffraction angle [rad] 

Fig. 7. Comparison of numerical calculated diffraction profiles 
(solid line) and measured intensities (dots), normalized to the 
primary beam intensity. The measured profile at U — 28.79 V 
(a) and U — 28.90 V (b) corresponds to the calculation with 
amplitudes Lp = 7.4° (a) and Lp = 12.6° (b). The corresponding 
deflections from measured /2//0 using (|37|) are 7.5° and 11.0°. 
The pattern wave length Adir = 48 pm, cell thickness 48.5 pm, 
t/c = 28.73 V. 



For measurements near the instability threshold, we 
use cell 2. We study the stationary director deflection am- 
plitude as a function of the control parameter e = /7//7c — 1 
by increasing the driving voltage gradually from a subcrit- 
ical value to a voltage above the threshold. The voltage is 
increased adiabatically slowly such that the director de- 
flection is always practically in equilibrium. Simultane- 
ously, the diffraction intensity at the second order reflex 
is recorded as a quantitative measure of the director field 
modulation. The constant scattering background is elimi- 
nated as above. The experiment is repeated with the op- 
posite direction of the field sweep, the diffraction intensity 
is recorded while the driving voltage is decreased with the 
same rate. The expected characteristics for a perfect for- 
ward pitchfork bifurcation \^ ip (x \fe. Considering relation 
(inZj), this would correspond to a quadratic dependence of 
the diffraction intensity at the second order reflex from 
the control parameter 



/2 cx (X e 



(38) 




"0.005 0.000 0.005 0.010 
Controll parameter e=^U/U-l 



0.015 



Fig. 8. Intensity at the second order peak within a sequence 
of increasing (o) and decreasing (A) driving voltage in cell 2. 
The fit to a function / oc (U — UcY yields a critical voltage Uc — 
28.84 V at sine wave excitation with v — 500 Hz. Due to a slow 
drift in conductivity the critical voltage is shifted to a higher 
value as in Fig.|3 In the insert, the same data are presented 
on logarithmic scale, an increase of fluctuation amplitudes in 
the subcritical voltage range is clearly observable. 



mically in the insert. The systematic deviations of data 
taken during up and down sweeps of the field, resp., are 
negligibly small. No hysteresis is found. The parameter 
Uc is fitted such to obtain best agreement with Eq. ((38j) . 
The measured intensities match the expected quadratic 
behaviour, and therefore we associate the fitted Uc with 
the threshold voltage. Close to the threshold and below 
[/c, the characteristics is covered by additional influences 
of noise. It clearly deviates from the prediction of Eq. 

indicated by the dashed curve. One of the possible 
reasons is that subcritical fluctuations of the modes close 
to the instability threshold piElET lead to an increased 
diffraction signal at the corresponding position. In addi- 
tion, the equations used in the hydrodynamic model use 
exact planar boundary conditions, while the cell actually 
has a small pretilt, typical for glass plates of sandwich 
cells with antiparallel rubbing. 



5.2 Measurement of growth and decay rates at 
periodic excitation 

The solution of the linearized differential equations Q 
yield an exponential growth or decay of ipt in an elec- 
tric field of constant amplitude. Equation (pTf)) connects 
this amplitude of the director deflection with the measur- 
able diffraction intensity, e.g. at the second order reflex. 
Both equations can be combined to 



/2(t)=/2(0)e^"^^* = /2(0)e^^- 



(39) 



FigurelHl shows the experimental results. For better visu- 
alization, data are presented in linear scale and logarith- 



where \i{E^v^kc) is the largest eigenvalue of the matrix 
product T+T~ and Aexp the experimentally determined 
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Time t [s] Time t [s] 



Fig. 9. Time evolution of the diffraction intensity at the sec- 
ond order reflex after changes of the field (20 Hz). The vertical 
axis shows the amplified voltage at the photo diode, which is 
proportional to the intensity I2 ■ The threshold [/c is 9 V in case 
of cell 1. In (a), the applied square wave voltage is switched at 
t = from zero to supercritical voltages. In (b), the excitation 
voltage is switched down to slightly subcritical values. 

growth rate from the intensity change at the second order 
reflex. The factor 4 considers the fourth order dependence 
of the scattering intensity from the director deflection am- 
plitude. Positive growth rates for Aexp > can be obtained 
in the experiment by recording the intensity change at the 
reflex after an electrical field E > Ec is turned on. For the 
measurement of Aexp < 0, the electric field is first switched 
to a supercritical value E > E^ where the convection pat- 
tern develops. Then, the field is suddenly changed to a 
value E < Ec and the intensity trace is recorded. Exam- 
ples for these procedures in cell 1 are shown in Fig.El In 
order to detect the fast changes we use a analog-digital- 
converter with a lower resolution (12 bit) but much higher 
sampling frequency of 1 kHz. The constant background is 
eliminated and the data are fitted to exponential func- 
tions in the middle of the detection range. Growth and 
decay rates of 200 measurements are depicted in Fig. QSl 
together with the eigenvalue Ai(£') calculated analytically 
from the material parameter in Table 1, and the factor 4 
from Eq. (|39|) . has been taken into account. There is good 
quantitative agreement with the linearized theory for elec- 
troconvection for driving voltages around the threshold. 



5.3 Trajectories under stochastic excitation 

If the deterministic voltage is replaced by a stochastic ex- 
citation sequence, the trajectory of the director deflection 
exhibits irregular changes in time. A statistical analysis of 
this phenomenon has been described in detail ^28729^ . The 
measurement of the time dependent diffraction intensity 
provides a convenient tool to study the trajectories of pat- 
tern amplitudes in real time. Figure ITTl shows the example 
of measured intensities at the second order reflex and the 
corresponding numerical simulation of the trajectory by 
solution of the differential equation (@|). 

The stochastic sequence in this experiment was a di- 
chotomous Markov process (DMP) with jump rate 160 s~^ . 




-0.2 -0.1 0.0 0.1 0.2 

Control parameter e = U/U^ -1 



Fig. 10. Growth and decay rates at different excitation am- 
plitudes (square wave). Growth rates (•) has been obtained 
from traces similar to those in Fig. and decay rates (A) 
from traces like those in Fig.l^J). The dashed line depicts the 
theoretical prediction 4iyXi{E) from a calculation of the largest 
eigenvalue of T^T~ in Eq. @. Material parameters are given 
in Table 1. 

Since the excitation sequence is synthesized with a com- 
puter, it is possible to use identical noise sequences in 
both experiment and simulation. In the bottom part of 
Fig. ^2 the realization of the stochastic driving process is 
shown. The numerical I{t) have been obtained from the 
(p{t) trajectories by use of (|37|) . Figure ITTl demonstrates 
that experiment and theory for stochastically excited EHC 
do not only agree on the statistical level when fundamen- 
tal scaling laws are compared, but even in details of the 
trajectories to a satisfactory degree, when we take into ac- 
count that the simulation cannot treat the involved addi- 
tive noise exactly but substitutes it by some average 28j. 
Of course, even repeated measurements of experimental 
trajectories are not exactly reproducible because of such 
additive (thermal) noise. Trajectories taken with the same 
noise sequence of the driving field differ in detail at small 
intensities, but above the noise level they are very similar 
and reproduce the simulated curve on average. 



6 Summary 

We have used laser diffraction as an detection method 
for direct quantitative determination of the amplitudes 
of the director field in nematic electroconvection. Instead 
of the evaluation of the complete diffraction pattern, it 
is sufficient to record the second order diffraction spot, 
which is mainly influenced by the phase grating generated 
by the director field. An analytic treatment of the ray 
propagation in the LC layer by Fermat's principle provides 
the qualitative relation /2//0 = const, x cp^ at the second 
order diffraction reflex for small pattern amplitudes. 

The proportionality constant can be derived from the 
material parameters using an analytical approach that 
considers at least the first non-linear term in light de- 
flection in the calculation of the optical path of individual 
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Fig. 11. Real time detection of director amplitudes in stochas- 
tically driven EHC, compared with a numerically simulated 
trajectory with the same realization of the driving process. The 
upper part (a) shows a measured intensity detected with the 
photodiode at the second order diffraction reflex (solid line), 
and a simulated curve (dotted line), both corresponding to the 
driving sequence depicted in (b). 



light rays passing the cell. A numerical calculation of light 
propagation, which does not use mathematical approxima- 
tions, except for the concept of ray optics, confirms the an- 
alytical result up to sufficiently large director deflections. 
The complete diffraction intensity profile calculated nu- 
merically is in good agreement with the profile measured 
experimentally. It has been shown that it is sufficient to 
use the simple relation Eq. ((57)) to determine the abso- 
lute value of the director deflection amplitude from the 
diffraction efficiency. We note, however, that a rigorous 
treatment of the problem of light propagation (for exam- 
ple by means of the FTDT method 44 ) is the only exact 
treatment of the optical problem. It has not been achieved 
yet. 

The derived quantitative relations have been used in 
three applications. In case of periodic sine wave excitation, 
the reported technique permits us to confirm the square 
root characteristics of the pattern amplitude in the pitch- 
fork bifurcation of the stripe pattern. The non perfect 
behaviour at subcritical values is attributed to thermal 
fluctuations and a slight sample pretilt. 

The agreement between analytically calculated growth 
and decay rates of the amplitudes of the director field 
and the measured light intensity at the second order re- 
flex of diffracted light shows that one has to be very cau- 
tious when growth/decay rates are determined from laser 
diffraction intensities. The decay of the optical signal goes 
with the 4th power of the director deflections, and conse- 
quently, time constants differ by a factor of four. 

The real-time quasi optical Fourier transformation 
of the pattern gives an easy access to fast changing mode 
spectra or amplitudes, e.g. in case of stochastic driving. 
We have demonstrated the direct correlation between the 
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Table 1. Material parameters in Eq. @ used in the calcula- 
tions. Experimentally data for Mischung 5 (last column) have 
been taken from jl4ll45) . measured conductivities correspond 
to the non-doped material. The unknown parameters and con- 
ductivities of the individual cells are obtained by fitting the 
corresponding threshold voltages and wave number character- 
istics for periodic ac driving to experimental data, see also 
Fig.ia 



driving electric field and the response of the director. Ob- 
servation of the whole diffraction image instead of the tra- 
jectory of one representative diffraction peak may provide 
access to the development of the mode spectrum and ac- 
cess to dynamic mode selection in the stochastic driven 
system. 

An important aspect in the experiment is the strong 
dependence of the diffraction image on small deviations 
from normal incidence. The problem of oblique incidence 
has been addressed first by Zenginoglou and Kosmopoulos 
PSj . In the appendix, we consider in detailed the depen- 
dence of the diffraction pattern from the angle of inci- 
dence of the laser beam. It is demonstrated that oblique 
incidence in general favors the reflexes of odd numbered 
order, which in first line reflect the amplitude grating pro- 
duced by the director field. 

The authors are particularly indebted to H. Schmiedel 
for helpful comments, discussions and critical reading of 
the manuscript. We acknowledge financial support from 
the Deutsche Forschungsgemeinschaft (Grant Be 1417/4 
and SFB 294). 
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A Diffraction at oblique incidence 

The experiments show that the diffraction profile has a 
strong dependence on a the tilt of the cell respective to 
the incident laser beam. A theoretical treatment with lin- 
earization has been published in j^. In case of oblique 
incidence, where the cell is tilted in the (x, z) plane, the 
symmetry ^ —0 is broken. The initial condition for Eq. 
(|TKj) is now r'{z) = tan 5^ where S is the entry angle of the 
laser beam into the LC-layer. We consider this in the nu- 
merical calculation of the beam propagation. In addition 
to this symmetry breaking, an extra phase AiIj{6) appears. 
Figure El illustrates the origin of the additional phase 
difference. For the primary beam at 6> = 0, Aip{0) van- 
ishes. The symmetry breaking leads to a slight shift of the 
diffraction spots, and more obviously, to a change of the 
relative intensities of even and odd order spots. The con- 
dition allowing for the additional phase for constructive 
interference on the 2nth even order spot (n > 0, > 0) 
is (see Fig. [T2t) 
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Fig. 13. Density plot from simulated diffraction profiles for 
oblique incidence and a cp = OA rad. The gray scale is log- 
arithmic in the intensities. A cell thickness of 25 pm and a 
pattern wave length of 33 pm assumed. The dashed line marks 
6= 
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where Aph is the wave length of the phase modulation of 
exiting light and the wave number of the laser light. 
The diffraction intensity calculated numerically for a di- 
rector modulation of = 0.4 rad as a function of Q and b 
is depicted in the density plot of Fig. ^1 The most obvi- 
ous result of the numerical calculation are the quantitative 
changes of the diffraction intensities with the cell rotation 
angle. Whereas the odd order maxima, which are mainly 
generated by the amplitude modulation, exhibit a mini- 
mum in the non-tilted cell and increased with a slight tilt 
of the cell, even order spots show qualitatively opposite be- 
haviour. In the experiment we record the complete diffrac- 
tion image on a diffusely scattering screen with a CCD 
camera (see Fig.|3| and scan the line y = (see Fig^t) 
from the digital image sequence. The camera gives only a 
qualitative picture of the intensities, not an exactly linear 
representation, but qualitative agreement with Fig. ^1 is 
clearly acknowledged. We remark that the smallest diffrac- 
tion angle On{S) for a given order n is not reached at (5 = 0, 
but 

Its position is depicted as dashed line in Fig. [T^ and coin- 
cides with the largest amplification/attenuation, resp., of 
the diffraction peak intensities. 



Fig. 14. Density plot of the measured diffraction profiles in 
inverse gray scale for a tilted cell. The profiles are taken in the 
midplane y — oi the diffraction images of a stable convec- 
tion pattern. The tilt of the cell to the incident beam leads 
to modulations of relative intensities and positions for all re- 
flexes. One acknowledges the symmetry {—S, —0) ^ {S,0). The 
dashed line marks ^ = t:^ as in Fig.ini 



